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Abstract. The mass distribution of the rich cluster of galaxies C10024-I-1654 has frequently been used to constrain 
the nature of dark matter yet a model consistent with all the observational data has been difficult to construct. 
In this paper we analyse the three-dimensional structure of this cluster using new spectroscopic information on 
~ 300 galaxies within a projected radius of 3/i~^Mpc. These data reveal an unusual foreground component of 
galaxies separated from the main cluster by 3000kms"\ We use numerical simulations to show that a high 
speed collision along the line of sight between C10024-I-1654 and a second cluster of slightly smaller mass can 
reproduce the observed peculiar redshift distribution. Such a collision would dramatically alter the internal mass 
distribution of the bound remnants, creating constant density cores from initially cuspy dark matter profiles and 
scattering galaxies to large projected radii, consistent with the observed distribution of galaxies in C10024-I-1654. 
The proposed scenario can reconcile the inferred mass profile from gravitational lensing with predictions from 
hierarchical structure formation models, while at the same time resolving the mass discrepancy that results from 
a comparison between lensing, velocity dispersion and X-ray studies. 
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1. Introduction 



of radius 106 ft. -'^ kpc that has been used to reconstruct 



The naLiiie uf the ubiquitous daik mat Lei leiuaiiis one uf 
the majoi m^isLuiies in Losiuolug^. The LiadlLiuiial vluw is 



the projected central m a ss distribution ( Kassiola et al 



that tile ddik maLLui is told, i. e. uuu-iuldlivisLiL at the 
decoupling epoch, and only subject to gravitational and 
weak interactions. Structure formation with this type of 
dark matter has been investigated extensively through nu- 
merical simulations over the past 25 years. High-resolution 
simulations of individual cold dark matter halos have 
shown that the mass distribution in CDM halos over a 
wide range of total masses, from dwarf galaxy to clus- 



1992|; ISmail et al. 199^ ; fPyson et al. 1998| ; [Broadhurst 



et al. 200C ). Constructing a mass model for this cluster 
that is consistent with all the observations has proved 
very difficult, however. The total projected mass enclosed 
within the arc radius is fairly well constrained and yields a 
value of about 1.6 x 10^^ h^^ ^^0- To reproduce this mass 
with an NFW type profile would require a cluster with 
characteristic velocity dispersion larger than 2000 km s""'^ 
Oroadhurst et al. 200d|; [Shapiro fc Ihev 2000|). This is 



ter scales, should follow a universal profile (Navarro et al 



1997, NFW), in paiticulai, the mass distiibutiun is found 
Ttsp with logarithmi c slope of "1.5 



to hav e a central 



much higher than observed in this or indeed in any other 
cluster of galaxies. The mass reconstruction by Tyson 
et al. (199S), based on the detailed structure of the grav- 



([Moore et al. 1998| ; [Ghigna et al. 2000[ ). Recent obser- 
vational results, however, have challenged this simple pic- 
ture. 

The rich cluster of galaxies C10024-I-1654 at .2 = 0.395 
features a spectacular quintuple gravitational arc system 
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itationally lensed background spiral galaxy, shows a flat 
core in the projected central mass distribution with a core 
radius of 35/i~^kpc. This apparent contradiction to the 
predictions of standard CDM simulations has prompted a 
number of authors to investigate alternatives to the clas- 
sic CDM scenario, e. g. warm dark matter, self-interacting 
([Spergel fc Steinhardt 2000| [Hogan fc 

2ooq ). 



Dalcanton 2000; Moore et al 
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As more detailed observations of individual clusters of 
galaxies are compiled, combinations of X-ray imaging and 
spectroscopy, velocity dispersion measurements and lens- 
ing mass maps frequently reveal that clusters that were 
thought to be simple relaxed objects are actually more 
complex systems, frequently undergoing mergers or gen- 
erally showing signs of substructure and deviations from 
dynamical equilibrium. Evidence for recent mergers or ac- 
cretion is currently accumulating through Chandra ob- 
servations in the form of merger shocks (e. g. Markevitch 
& Viklhnin 20011; iMarkcvitch ct al. 2001[) or cold fronts 



-4000 



V [km/s] 

-2000 2000 4000 



(Markevitch et al. 200C; Vikhlinin et al. 2001; Mazzotta 



et al. zOOl). Mergers and substructure have also been in- 



voked to explain the discrepancy between mass estimates 
from different methods that is observed in many clusters 
(iMiralda-Escude fc Babul 1995|; |Wu et al. 199^). 



C10024-I-1654 is an example of a cluster where the high 
mass inferred from the strong lensing mass reconstruc- 
tions is at variance with the fairly low X-ray luminosity 
and temperature, that indicate a total mass for the clus- 
ter which is a factor 2 to 3 smaller than the lensing mass 
( ^oucail et al. 2000| ). Both the galaxy distribution on the 
sky and the X-ray morphology are regular and by them- 
selves compatible with a relaxed massive cluster, a notion 
which was further supported by the high galaxy velocity 
di spersion of 1200 km s~ ^ fo und in the rcdshift su rveys 

(19991, 



here- 



of Dressier fc Gunn (1992 ) and Dressier et al. 
after D99). 

Paper I) we presented a new 



In Czoske et al. (2001 



catalogue of 650 measured redshifts (including those from 
D99) in a wide field around C10024-hl654. 300 galaxies in 
the catalogue have redshifts in the range 0.37 < z < 0.42, 
i. e. are cluster members or lie in the immediate neigh- 
bourhood of the cluster, and the redshift distribution of 
these galaxies was found in Paper I to be distinctly bi- 
modal. In the present paper we analyse the distribution 
of the galaxies in redshift/real space and present a merger 
scenario which is able to reconcile all the observations. 
Section || recalls the data used and describes the distribu- 
tion of the redshifts globally and as a function of position 
within the survey field, as well as the distribution of spec- 
tral types based on equivalent widths of important lines. 
In Section |[ we develop the merger scenario and use nu- 
merical simulations to show that it is possible to explain 
the peculiar three-dimensional distribution of the galaxies 
in C10024-I-1654. In Section ^ we show how the observed 
structure is able to explain and resolve the discrepancy 
of mass estimates of C10024-I-1654 derived from different 
observations. Finally, in Section ^ we summarize and con- 
clude. 

Throughout this paper we use a Hubble constant Hq = 
100 ft,"""^ kms""'^ Mpc""'^ and assume an Einstein-de Sitter 
Universe with J^m ~ 1 and n\ = 0. At the redshift of 
C10024+1654, z = 0.395, 1" corresponds to 3.195 /i^^ kpc. 



2. The cluster environment 
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Fig. 1. Redshift histogram for 300 objects in the neigh- 
bourhood of C10024+1654 (0.37 < z < 0.42). The attribu- 
tion of galaxies to components A (dark grey) and B (light 
grey) was done by inspection of Fig. ||. Whereas the lower 
and upper redshift limits of components B and A respec- 
tively are fairly secure, the exact demarcation between 
the two components is arbitrary to some extent and was 
chosen at Aw — — 1500kms~^ with respect to the central 
redshift of component A, Eq. (0). 



2.1. The data 

In Paper I we presented a new catalogue of photometric 
and spectroscopic data for 679 objects in a wide field of 
21x25arcmin^ around C10024-t-1654. At the redshift of 
the cluster, z = 0.395, this corresponds to 4 x 4.8 h^^ Mpc^ 
(the true survey field is delimited by an irregular polygon 
inscribed in this rectangle and covers 16.8 ft-^^ Mpc^, see 
Fig. |[ or Fig. 7 in Paper I). The catalogue lists equatorial 
position, redshift, V magnitude, V — I colour, equivalent 
widths for [O ii] A3727, [O ill] A5007, Ha (where within the 
wavelength range), H/3 and US, as well as the strength of 
the 4000A break. Here we are only concerned with the 
650 objects having measured redshifts, of which 581 were 
observed by u s at CFHT and WH T, with the remaining 
69 taken from Dressier et al. (1999| ). See Paper I for more 
technical details on the data used in this paper. 



2.2. Redshift distribution 

Fig. |l| shows a histogram of the 300 redshifts that lie in 
the range 0.37 <z< 0.42. The redshift distribution of the 
cluster galaxies is clearly bimodal, showing two peaks at 
z — 0.381 and z — 0.395 respectively; these peaks contain 
283 galaxies. 

The larger peak at z = 0.395 (0.387 <z< 0.402, here- 
after referred to as component A) contains 237 galaxies, 
is fairly regular and resembles a Gaussian distribution as 
expected for a relaxed, virialized cluster of galaxies. The 
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smaller foreground peak at z = 0.381 (0.374 <z< 0.387, 
component B), by contrast, seems too wide for the small 
number of 46 redshifts contained in it. This impression 
is confirmed by Fig. |^, which plots the redshift for each 
galaxy versus its projected distance from the centre of 
C10024+1654|]. The distribution of the galaxies in the 
main peak is symmetrical with respect to the central 
redshift line, whereas the distribution of the foreground 
galaxies is roughly constant at a rest frame velocity of 
~— 3000kms~^ at radii larger than 3' (600 kpc), but 
turns off towards smaller relative velocities to merge with 
the main distribution at smaller projected distances. 

The 17 galaxies in the peak at z = 0.407 are more 
widely dispersed across the survey field and although we 
cannot rule out a connection with C10024+1654, it seems 
more likely, in the light of the scenario developed in Sect. 
^, that they are part of the surrounding field galaxy pop- 
ulation. 

It is remarkable that we can trace the main cluster 
as well as the foreground structure out to the edge of 
the survey field (there remains only one object in the 
catalogue outside the range depicted in Fig. ||, a back- 
ground galaxy). The virial radius for C10024+1654 should 
be around l/i^^Mpc ( pirardi fc Mezzetti 200 1| ), which 
means that we detect a coherent structure out to three 
times the virial radius. 

Redshift distribution inside 200" 

From Fig. || and Fig. ^ (upper panel) it is obvious that 
the central region of C10024-t-1654 is highly perturbed and 
it is impossible to separate components A and B within 
~ 200" (640 Mpc) from the projected cluster centre. 
However, the redshift distribution for the central galaxies 
is strongly skewed towards negative velocities: except for 
two galaxies that are fairly isolated in redshift space, the 
redshift distribution at positive velocities is effectively cut 
off at ~ 1500 km s^^, whereas the distribution at negative 
velocities extends to beyond — 2500 km s^^, branching 
off into component B at ~ 500 h^^ kpc from the projected 
cluster centre. Quantitatively the skewness of the distri- 
bution of 161 galaxies within 200" (0. 37<z<0.42) is 0.5 4, 
where we use the definition given by Press et al. (1992): 



central redshift distribution is already apparent in the his- 
togram given by Dressier et al. (199^ ) (see also Section ||). 



Skew(a;i 



■ xn) 



1 ^ 



(1) 



The probability that a sample of size A^ = 161 drawn from 
a normal distribution shows a skewness larger than that 
measured in C10024-t-1654 is about 0.3%. The Shapiro- 



Wilk normality test ( Shapiro fc Wilk 1965 ) rejects the hy- 
pothesis that these data points are drawn from a Gaussian 
parent population at ^ 99% confidence. The skew in the 

^ The coordinates in the catalogue are given relative to 
02000 = 00*^26^35': 70, ^2000 = 17°09'43"06. Originally, this 
was the position of galaxy 373, but a more accurate astromet- 
ric analysis shifted the reference point by 0'.'79 to the south-east 
of this galaxy. 



Redshift distribution outside 200" 

In the distribution of the galaxies outside 200" (Fig. || 
and Fig. ||, bottom panel) the two modes are clearly sep- 
arated. The mean redshift of the 71 galaxies belonging to 
component A at projected distances 200"<i?<500" is 

il = 0.3946 ±0.0007 . (2) 

In Fig. H we show the two velocity dispersion profiles for 
galaxies in component A with either positive or negative 
velocities with respect to the mean redshift za, calcu- 
lated in a sliding bin containing 30 galaxies. Outside 3' 
(600 kpc) both profiles are fiat at the same level of 
(J ~ 600 km s^^. The profile for galaxies with negative ve- 
locities (i. e. those moving towards us, Fig. ^ rises all 
the way to about 45" (140 kpc), reaching a velocity 
dispersion of nearly 900kms~^ in the innermost bin. The 
dispersion profile of the galaxies with positive velocities, 
by contrast, drops back to a value of about 600 km s^^ af- 
ter having reached a maximum value of about 800 km s^^ 
at 2' from the centre. Although the galaxy velocities used 
in Fig. H have been restricted to |u| < 1500 kms~^, the 
different behaviour of the profiles for galaxies with posi- 
tive and negative velocities respectively reflects the skew 
in the central redshift distribution and indicates the pos- 
sible presence of a bulk motion component towards us in 
the central galaxy distribution. 

The exact value of the velocity dispersion of the 71 
galaxies within the same distance range as used for Eq. 
(||) and velocity |u| < 1500 km s"' 
redshift za is 



relative to the central 



561^^^ kms" 



(3) 



Bot h ZA and a a were c omputed using the biweight estima- 
tor ( Beers et al. 1990 ), the errors (95% confidence level) 
were estimated by bootstrap resampling. 

The skewness of the distribution of these 71 redshifts 
is —0.12, which should be compared to the standard de- 
viation of skewness for samples of this size drawn from a 
normal distribution, 0.27. The Shapiro- Wilk test does not 
provide evidence for deviations from a Gaussian distribu- 
tion. 

The spatial positions of the galaxies in the foreground 
component B show no preference for any direction (Fig. 
^, middle panel). The formal velocity dispersion for com- 
ponent B in the same radial distance range as used in Eq. 
(||) is (Tb 554^304^"^^"^ based on 15 redshifts. Given 
that the number of galaxies in component A is almost five 
times as large as the number in component B it is hardly 
conceivable that both components should have the same 
velocity dispersion - their mass-to-light ratios would be 
extremely different. Component B is certainly not a viri- 
alized group or cluster. 

From the distribution of the external galaxies alone 
one might still suspect the presence of a loose aggregate 



4 



O. Czoske et al.: Spectroscopic survey of C10024+1654: II. 



of galaxies physically disconnected from C10024+1654, but 
well-aligned with the line of sight; however, in the centre 
there are no galaxies which would correspond to the z = 
0.38 peak visible in the external distribution (Fig. ||). If 
the foreground component B were an independent loose 
system, it would show a deficiency of galaxies right in front 
of the centre of C10024-I-1654, which seems unlikely in the 
absence of interaction between the systems. This leads us 
to suspect that the two components are in fact physically 
connected and that the blue tail of the skewed central 
distribution originates from the same physical system as 
the foreground component B. 

Fig. 1^ shows the theoretical escape velocities for masses 
of 1, 3 and 5 X lO^M© enclosed in a sphere of radius R. 
Although these lines are merely indicative due to the as- 
sumption of spherical symmetry and the unknown tangen- 
tial velocity components, it is safe to say that the galaxies 
in component B outside, say, 600/i~^kpc are not bound 
to the main cluster component, unlike the negative tail of 
the skewed central galaxy distribution. 

The galaxy density map of the main component A 
(Fig. 1^) shows an extension of the galaxy distribution to- 
wards the northwest to a distance of 100" — 150". This ex- 
tension is visible in images of C10024-I-1654 and we show 
here that the galaxies contained in it lie at the redshift 
of the main cluster. C10024-I-1654 might therefore well be 
comprised of three distinct subsystems: the main cluster, 
component A; a foreground cluster or group, component 
B, which is interacting with the main cluster; and a third 
group or small cluster, which is falling onto the main clus- 
ter in the plane of the sky. 



2.3. Distribution of spectral types 

In Tables |^ and |^ we classify the galaxies into spectral 
types, defined as in [Balogh ct al. (1999 ), using our mea- 
sured values for the equivalent widths of [O ii] , H (5 and the 
4000A break (see Paper I for details on how these values 
are defined and measured). 

We consider four samples of galaxies according to red- 
shift and projected spatial position; the sample bound- 
aries for samples "Centre", "A" and "B" are marked by 
the dotted lines in Fig. ^. The boundary between "Centre" 
and "B" is to some extent arbitrary; replacing the slanted 
line by a simple velocity cut at v — — 1500kms~^ does 
not significantly change the numbers in Tables [| and ||, 
however. "Field" includes all the galaxies in the sample 
with redshifts < z < 0.55 without the cluster galaxies, 
0.372 < z < 0.402 (the upper limit here is the same as for 
the field sample in Balogh ct al. (1999)). 

Both the classifications by [On] vs. US and by KS 
vs. the 4000 A break show an excess of blue star- 



seems to provide support for the hypothesis that compo- 
nent B is a loose galaxy overdensity with spectral charac- 
teristics typical for field galaxies. Alternatively it is possi- 
ble that the star-formation activity in the galaxies belong- 
ing to component B has recently been (re) activated and 
enhanced through interaction with a denser environment. 

We include in Tables |l| and |^ corresponding numbers 
from Balogh et al. (1999| ). These numbers are not strictly 



forming galaxies (star-forming/short starburst, SF/SSB, 
and emission- line galaxies, A-|-em) in component "B" as 
compared to the cluster centre and possibly even over the 
field. The distribution of spectral types in component "A" 
is between a typical cluster and field populations. This 



comparable, because the selection criteria used for the cre- 
ation of our catalogue are quite different from those used 
by Balogh et al. (1999[ )p| this explains the discrepancies in 
particular in the US vs. the 4000 A break classification, 
although the general trends are the same. The agreement 
in the [O il] vs. H S classification on the other hand is very 
good, even quantitatively. 

3. A high speed collision? 

In general, a redshift difference between well-separated 
clusters at cosmological distances is a combination of a 
cosmological redshift difference due to the proper distance 
between the clusters along the line of sight and a Doppler 
shift due to the relative velocities of the two clusters. In 
general, these two effects cannot be disentangled uniquely, 
and there might be several ways to reproduce a given red- 
shift distribution. In the following we develop a scenario 
that can reproduce the observed redshift distribution in 
C10024-f 1654, in particular the skewed central distribu- 
tion, the presence of two well-separated modes away from 
the projected cluster centre extending to at least three 
times the virial radius. 

Assume that a group or cluster has undergone a near 
radial collision with C10024+1654 and we arc now viewing 
the merging system along the direction of the impact. In 
fact, several characteristic parameters of a collision sce- 
nario are fairly well constrained by the observed redshift 
distribution. The relative velocity of the two subclusters 
is given by the redshift difference of components A and 
B and is thus of order 3000 km s~^. Those parts of the 
smaller cluster B that pass through the centre of A are 
more strongly decelerated in the direction of the collision 
than its outer parts due to stronger dynamical friction 
and they therefore show a smaller redshift difference. A 
tidal gravitational shock during the crossing scatters the 
outer galaxies of both clusters to large projected distance, 
thus increasing the extent of the cluster halos beyond their 
nominal virial radius. If the transverse velocity imparted 
on a galaxy that is now found at a radial distance of 
3h~^ Mpc is 1000 km s~^, then the time since core crossing 
is about 3 Gyrs, and the separation between the cores of 
the two clusters is roughly 5 Mpc. 

The main unknown in this scenario is the mass ratio 
between the subclusters. There are two possible versions 
of this collision: 



^ [Balogh et al. (1999| ) use data from the CNOCl cluster red- 
shift survey, which has limiting Gunn r magnitudes between 
20.5 and 22.0, de pending on the redshift of the observed cluster 
(|Yee et al. l"996[). 
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Fig. 4. Galaxy number density maps for (a) component A (0.388 < z < 0.402, 237 galaxies), (b) component B (0.37 < 
z < 0.388, 46 galaxies) and (c) the "field" (0.1 < z < 0.35, 141 galaxies), estimated with the generalized nearest 
neighbour method (Silverman 1986) with 10 neighbours for map (a) and 5 neighbours for maps (b) and (c). The maps 
are divided by the completeness map (see Paper I) and smoothed with a Gaussian with a — 30". The grey scales are 
the same for all three maps. The lowest density contour lines are at 2% and 5%, normalized to the maximum density 
in the cluster component A (left panel). The remaining contours are spaced in steps of 10%. For densities below 50% 
the contours are drawn in black, for higher densities in white. The 50% contour is marked by a bold white line. 



Table 1. Distribution of spectral types: classification by [On] 
from 



H(5 (the definition of the spectral types is taken 



Balogh ct al. (1999) and given in column 2). The fraction of galaxies of a given spectral type is listed as the 



percentage of the total number of galaxies in the subsample (the total number is given in parentheses). See text and 
Fig. H for the definitions of the samples. The numbers are the percentages p of galaxies belonging to each class in each 
sample. The errors are the standard deviations Np{l — p) of the normal distribution asymptotic to the binomial 
distribution with sample size N and probability p for large N. 













This work 




Balogh et al. (1998) 


Spectral Type 


Definition 




Centre (183) 


A (64) 


B (35) 


Field (275) 


Cluster 


Field 




H5< 


-5 A, [On 


<5A 


2±1 


3±2 





4±1 


4 


2 


A+em'' 


H5< 


-5 A, [On 


>5A 


3±1 


6±3 


11±5 


3±1 


4 


10 


SF/SSB'^ 
passive"^ 


H(S> 


-5 A, [On 


>5A 


23±3 


42±6 


57±8 


53±3 


24 


47 


H(5> 


-5 A, [On 


<5A 


72±3 


48±6 


31±7 


40±3 


67 


41 



''spectrum dominated by K and A stars; ''spectrum dominated by A stars, with emission lines; ''star- forming/short starburst; 
"^red, passively evolving galaxies 

Table 2. Distribution of spectral types: classification by US and the 4000A break. The structure of the table is the 
same as for Table |l|. 











This work 




Balogh et al. (1999) 


Spectral Type 


Definition 


Centre (148) 


A (60) 


B (33) 


Field (225) 


Cluster 


Field 


bHDS'' 


H<5< 


-5A, br<1.5 


3±1 


3±2 


9±5 


7±1 


5 


7 


rHDS'' 


H<5< 


-3A, br>1.5 


8±2 


17±5 


9±5 


5±1 


10 


11 


SF/SSB" 
passive'^ 


H(5> 


-5A, br<1.5 


34±3 


50±6 


67±8 


56±3 


12 


33 


H5> 


-3A, br>1.5 


55±4 


30±6 


15±6 


32±3 


74 


45 



''blue H5 strong; ''red H5 strong; ''star-forming/short starburst; ''red passively evolving galaxies 



Fig. 6. The initial and final particle configurations used to simulate the two colliding clusters. The grey scale indicates 
the local density of dark matter within the box of length 12 Mpc. 
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(i) A small group has passed through the core of 
C10024+1654 and has been completely disrupted and 
scattered by the cluster's tidal force. In this case the 
main body of C10024+1654 is unperturbed and we do 
not observe a concentration of foreground galaxies be- 
longing to the original group because the system has 
been completely unbound and is scattered to large dis- 
tances perpendicular to the main cluster. 

(ii) A massive cluster of approximately 50% of the mass 
of C10024-I-1654 collided and passed through the core 
of the main cluster roughly 3 Gyr ago. This collision is 
insufficient to completely disrupt the impacting cluster 
but its outer galaxies are scattered to large projected 
distances. The remaining bound core of the colliding 
cluster undergoes sufficient dynamical friction to turn 
its orbit around and in redshift space it appears to lie 
at a similar distance to C10024-I-1654. 

Both of these scenarios can reproduce the main features 
of the observed redshift distribution. However, the mas- 
sive collision scenario is more interesting because of the 
possibility of resolving both the mass discrepancy prob- 
lem and the conflict between the central density structure 
of C10024-t-1654 and predictions of hierarchical clustering 
models. 

The mass 
al 



et 



reconstructions of C10 024 

1998| ; [Broadhurst et al. 2000|) using 



1654 ([Tyson 
the posi- 



tions of multiple images of a gravitationally Icnsed back- 
ground galaxy show a projected mass of 1.3 x 10^^ h^^ Mq 
within 106/i~^kpc and a central surface mass density of 
7900 HMq pc^^. For a singular isothermal halo extending 
to ?'2oo/kpc ~ \/2 f id/(kms^^) this implies a characteris- 
tic velocity dispersion larger than 1500 km s~^ in order to 
explain the projected mass. 



As Shapiro fc Iliev (2000 ) point out, the situation is 



worse in the context of hierarchical clustering models. 
Dark matter halos in cold dark matter (CDM) type mod- 
els are shallower than isothermal in their centres. In order 
to reproduce the observed projected mass a CDM halo 
with velocity dispersion of 2200 km s~^ is required. This is 
inconsistent with the observations presented in this paper. 



A second problem for hierarchical models is that Tyson 



et al. ( L998) infer a very shallow central density profile for 
this cluster, much flatter than the cuspy density profiles 



found for clusters in CDM type models ( Ghigna et al 
2000()~|The central core of C10024+1654 has frequently 
been used to constrain the nature of dark matter and 



to arg ue for alternative candidates to CDM ( Spergel fc 
Stei nhlrdt 2GO0| ; |Hogan fc Dalcanton 2000| ; [Moore et a~ 
20000 

A high speed encounter between two similar mass clus- 
ters can explain all of the observations presented in Sect. 
^ and reconcile the discrepancy between mass estimates 



derivec for C10024+1654. We now explore this scenario 



using high resolution numerical simulations of colliding 
dark matter halos to study the evolution of the mass dis- 
tribution. We construct two equilibrium CDM halos with 
virial masses 9.5 x 10^^ Mq and 5.0 x IO^^Mq with peak 



circular velocities Wpoak = 1600 and 1300 km s~^ and con- 
centrations c — 5 and c = 7 respectively. Their initial sepa- 
ration is 3Mpc and relative velocity is — 3000kms~^. The 
particle mass is set to 5 X 10^ Mq and we use an equivalent 
Plummer force softening of 5 kpc. 

After the collision we find that the outer regions of 
the smaller cluster have become unbound and are stream- 
ing radially away from the impact location. The impulse 
velocity perpendicular to the encounter is of the order 
1000 km s""'^. Snapshots of the initial and final times are 
shown in Fig. ^. The peak circular velocities of the bound 
components have fallen to 1430 and 1140kms^^ whilst 
the central ID velocity dispersions have reduced to 930 
and 710 km s^^ from initial values of 1033 and 883kms~^ 
respectively. At the final time we find a total mass within 
a cylinder of radius 106 kpc of 8.0 x IO^^Mq. 

The redshift distribution of 10 000 randomly selected 
particles at the final time is shown in the right hand panel 
in Fig. 1^. This is not a perfect match to the observational 
data but the main features are present: we see a fore- 
ground component of "galaxies" that span large projected 
distances from the central region of C10024-I-1654. This 
component is separated by « 3000 km s~^ in redshift space 
from the main component. The smaller cluster bound core 
is moving away from the main cluster at « 1000 km s~^ 
- the reduction in speed from the initial velocity is due 
to dynamical friction. A better agreement would result if 
the impacting cluster was slightly more massive in which 
case it would suffer more friction and would be seen nearly 
at rest compared to the main component. Also note that 
we are plotting "dark matter particles" not "galaxies". 
Of course the initial configuration shown in the left panel 
in Fig. also displays a bimodal distribution of redshifts. 
However, this configuration is more symmetric than either 
the final simulated configuration or the observed redshift 
distribution; also, the radial extent, in particular of the 
foreground (in redshift space) component, is much smaller 
than observed. 

The energy from the impulsive tidal shock has been 
transferred into kinetic energy of the particles, heating 
and expanding the cold central cores of the clusters. This 
leads to a flattening of the physical and projected den- 
sity profiles. In Fig. ^ we show the projected surface mass 
density profile before and after the encounter. The initial 
profiles are cuspy, CDM type density profiles whereas the 
final profiles have nearly constant density cores in good 
agreement with that inferred from the mass reconstruc- 
tion ( Tyson et al. 1998 , their Figure 4). The central surface 
mass density is 4000 M© pc~^ which is very close to the 
value obtained by Tyson et al. for Hq — 50kms~^Mpc~^. 
The total projected mass within the central 106 kpc is 



about 30% lower than obtained by Broadhurst ct al 
(2000D which could be reconciled by using a more massive 



encounter and/or properly including the baryonic matter. 
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4. Comparison to other observations 



4.2. X-ray observations 



The galaxy distribution of C10024+1654 in redshift space 
and on the sky, as described in Sect. ||, provides strong 
hints that this is not a simple relaxed system. How do 
the skewed central velocity dispersion and the presence of 
a foreground component in redshift space affect the in- 
terpretation of other observations and in particular mass 
estima tes of this system? 



The X-ray observati o ns from Rosat/HRI and ASCA 



( ISoucail et al. 2000 



Bohring em- 
selves show no indication of anything other than a fairly 
small cluster of galaxies. Given the measured values 

keV and Lx = 6.8 x 



5.7: 



for C10024-H654 Tx 

10^3 erg s"^ ( ISoucail et al. 2000D , the g as does not 



-4.9 



seem too far away from the Lx — Tx relation of Markevitch 



1998 ); it is slightly too hot for its luminosity. The mor- 



4.1. Galaxy dynamics 

Dynamical mass estimates rely on the assumption that 
the galaxies are in virial equilibrium in the cluster's gravi- 
tational potential well and that their velocities are purely 
random. 

The large velocity dispersion of a ~ 1200 km s^^ found 
by [Dressier fc Gunn (199"2|) and [Dressier et al. (1999D 



(hereafter D99) was based on galaxies within 5' from the 
projected cluster centre. As shown in Fig. ||, there is no 
direct evidence for the bimodality of C10024-I-1654 when 
attention is restricted to the central regions of the cluster. 
This is still true with the larger number of galaxy spec- 
tra used here: Whereas our catalogue contains about 85% 
more r cdohifto in thio area than D99, the general ohapo of 



the hi stogram io the oamc. Both hiotogramo arc otrongly 
skewed towards l ow redshifts, thus ind icating the presence 
of substructure ( Ashman et al. 1994 ) or a bulk velocity 
component. Note that the bimodality of the redshift dis- 
tribution only becomes apparent when galaxies at larger 
projected distances from the cluster centre are included. 
In this case it is therefore not the number of galaxy red- 
shifts but their distribution over a wide field that provides 
the clues to the complexity of this system. 

The formal velocity dispersion for the 193 galaxies 
within 5' (Fig. ^)is cTccnt = 1050kms~^. If there is indeed 
a bulk velocity component present in the central velocity 
distribution, then mass estimates based on the formal cen- 
tral velocity dispersion will overestimate the true cluster 
mass. 

A better dynamical mass estimate might be obtained 
from the velocity dispersion at larger radii, where the two 
components A and B are clearly separated. At projected 
distances > 3' the velocity distribution of component A 
is regular with a dispersion of ^ 600 km s^^. If we inter- 
preted this value as the velocity dispersion of a relaxed 
cluster, it would have only about a quarter of the mass 



phology is regular and the surface brightness profile well 
fitted by a beta profile with a surprisingly small core ra- 
dius of 33/i~^kpc. Whereas no significant substructure is 
required to model the X-ray observations of G10024-I-1654, 
they do not provide evidence against the collision scenario 
either. 

The X-ray emission is difficult to predict in the con- 
text of two clusters projected along the line-of-sight sev- 
eral Gyr after a head-on collision. During the collision, 
hydrostatic equilibrium in the gas component of the clus- 
ter (s) breaks down and the X-ray luminosity and emission- 
weighted temperature fiuctuate considerably (by up to 
a factor of 10 for Lx) and rapidly as a consequence of 
the formation of shock waves a nd repeated exp a nsion and 
compression of the core gas (Takizawa 1999 ; Ricker fc 
Sarazin 200l[ ; [Ritchie fc Thomas 200 1[ ). After several Gyr 



however, shocks have dissipated and the gas settles down 
to an equilibrium configuration. The cited hydrodynamic 
simulations of cluster mergers consider collision speeds of 
about 1000 km s~^; in this case the dark matter cores sep- 
arate to distances of only a couple of Mpc before turn- 
ing around and eventually merging into one clump. In a 
high-speed collision, however, the coUisionless dark mat- 
ter cores of the two clusters separate to large distance 
before turning around. The gas, due to its coUisional na- 
ture should experience stronger interaction during the first 
crossing and thus behave quite differently from the dark 
matter/galaxies and possibly also from what low-speed 
collision simulations predict. Ricker fc Sarazin 1200% ob- 



previously estimated for C10G24-I-1654 (Schneider et al 
1986| ), below mass estimates from gravitational lensing 



(which are roughly consistent with a high velocity disper- 
sion) and even below mass estimates from X-ray obser- 
vations. In the context of the collision scenario presented 
in Section ^ on the other hand, the galaxies at large pro- 
jected distance from the cluster centre are also affected by 
the collision and cannot be used to derive a mass estimate 
based on the assumption of dynamical equilibrium. 



serve a slight segregation between the gas and the dark 
matter cores in their simulations. This effect should be 
more pronounced in a high-speed collision. 

If we accept the collision scenario then a mass estimate 
from the X-ray observations is not possible. Detailed in- 
clusion of a hydrodynamic treatment of the gas compo- 
nent in simulations of high-speed collisions is necessary to 
understand the behaviour of the X-ray emission in this 
case. High resolution X-ray imaging with Chandra has 
recently been obtained (P. I. Hattori) but has not been 
published yet. 

4.3. Gravitational lensing 

Gravitational lensing by clusters of galaxies does not 
rely on the matter in the clusters being in dynamical 
equilibrium. However, gravitational lensing measures the 
weighted integral of all the mass between the observer and 
the source, and the interpretation of the measured mass 
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value thus depends on the detailed distribution of mass 
along the line-of-sight.The skew in the central redshift dis- 
tribution in C10024-I-1654 is an indication of bulk motion, 
possibly due to a high-speed collision between two fairly 
massive clusters of a mass ratio of roughly 2:1, or by the 
passage of a small group through the core of C10024-I-1654. 
In the latter case the mass determined by strong lensing 
models would indeed be indicative of the mass of the main 
cluster; however, in this case one would expect the X-ray 
emitting gas to be largely unperturbed by the collision, 
and one would be left with the discrepancy between the 
mass estimates derived from lensing and X-rays. In the 
former case, the lensing mass would be the sum of the 
two cluster cores and therefore too large if interpreted as 
representing the mass of a single cluster. 

The projected mass profile for the post-merger system 
seen in our simulations shows two characteristics which are 
testable by gravitational lensing analyses of C10024-I-1654: 
At small scales of several lO's kpc, the profile shows a flat 
core which is consistent with strong lensing mass models 
( Tyson et al. 19981) , as shown in Sect. |[ On large scales 
the profile is significantly flatter than the initial CDM pro- 
file, falling only as r^^-^ out to ^3Mpc. Previous weak 
lensing analyses on C10024+1654 have not attempted a 
detailed reconstruction of the mass profile. Bonnet et al 
(1994 |)land |van Waerbcke et al. (1997[) derived shear maps 



a numerical simulation we show that it is possible to ex- 
plain the observed redshift distribution with a high-speed 
collision of two rather massive clusters of galaxies with a 
mass ratio of about 2:1. 

Apart from reproducing the spatial/redshift distribu- 
tion of the cluster galaxies, this scenario also produces a 
projected mass distribution which is very close to that de- 
rived by Tyson et al. (1998) from the quintuple arc system 
observed in C10024-I-1654. Note that the scenario was not 
designed to reproduce this mass profile. Since it is thus 
possible to produce a mass distribution with a flat core 
from the merger of two CDM halos this eliminates one of 
the main arguments against simple non-interacting cold 
dark matter as the dynamically dominant com ponent in 
clusters of galaxies ( ^pergel fc Steinhardt 200C| ). 

The distribution of spectral types in the foreground 
component B is more akin to the general field population 
than to a cluster population. In the context of our colli- 
sion scenario, these would correspond to the outer regions 
of the smaller cluster which have become unbound dur- 
ing the impact. Even initially, these galaxies would prob- 
ably not correspond to a fully transformed cluster galaxy 
population. In addition, the impact might have triggered 
additional star formation in these galaxies. Remembering 
that C10024+1654 was one of the most distant clusters in 



from g: 'ound-based data on the north-east quadrant of the 
cluster, ISmail et al. (1997] ) were restricted to the small 



fiel d of view of HST/W FPC2. The shear pattern found 
by Bonnet et al. (199^ is compatible with an isothermal 
sphere profile, p cx r~'^ out to 3/i~^Mpc. 

A weak lensing analysis of newly obtained 
CFHT/CFH12k images in BVRI is currently under- 
way (Czoske et al. 2002, in preparation). Also, a sparsely 
sampled HST/WFPC2 mosaic consisting of 38 pointings 
out to a distance of 2.5/i^^Mpc from the cluster centre 
has recently been obtained; these data will also be used 



for a weak lensing analysis (Treu et al. 2001). These data 



which the Butcher-Oemlcr effect was detected (Butcher 
& Oemler 1984; Dressier et al. 1985), this may add a 
new view on this effect: The fact that ^ 40% of the 
bright galaxies in the cluster are emission line galaxies 
( [O ii] > 5 A) has a natural interpretation in the collision 
scenario. The spectral distribution of the galaxies in com- 
ponent A (and centre) are also perturbed by the collision 
with a possible excess of rHDS and SF/SBB (Tables |l]and 
H); rHDS galaxies may represent the result of a starburst 
induced during the early stages of the interaction. 

The relative velocity of the clusters of ~ 3000 km s^^ 
implied by the redshift distribution in C10024+1654 is 
very high. Observations and simulations find mean pe- 



will allow to probe the cluster's projected radial mass 



accurately to the edge of the survey field. 


& Oh 1996 


; Giovanelli et al. 1998 


IColberg et al. 2000 


). Still, coUiding 



Gibbons et al. 2001 



5. Discussion and conclusions 

In the present paper we have analyzed the galaxy distri- 
bution in the cluster of galaxies C10024-f 1654, based on 
~300 galaxy redshifts and projected positions. The clus- 
ter, which was previously regarded as a prototype of a 
massive relaxed cluster at intermediate redshift, turns out 
to have a fairly complicated structure, showing a strongly 
skewed redshift distribution in its central parts and two 
well-separated components at larger projected distances 
out to -3/i~iMpc. 

We argue that the blue tail of the central veloc- 
ity distribution and the foreground component originate 
from the same physical system and interpret the peculiar 
redshift-space distribution as the result of a high-speed 
head-on collision of two clusters of galaxies, the merger 
axis being very nearly parallel to the line-of-sight. Using 



IMpc ( ISarazin 200 1[ ). [Markevitch et al. (20011 ) have re- 



cently found a bow shock in the galaxy cluster 1E0657-56 
which implies a relative speed of 3000 to 4000 km s~^ for 
the collision. C10024-t-1654 is still an extraordinary but 
not an impossible system. 
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Fig. 3. Redshift histograms for objects inside (upper 
panel) and outside (bottom panel) a radius of 200" (corre- 
sponding to 640 kpc) from the projected cluster cen- 
tre. In the central redshift distribution there is no corre- 
spondence to the distinct peak at z = 0.38 visible in the 
external distribution. 



Fig. 2. Redshift z plotted against angular distance R from 
the projected cluster centre for the galaxies around the 
cluster redshift. The left axis expresses redshift as rela- 
tive velocity with respect to the mean redshift of com- 
ponent A (Eq. the right axis as proper line-of-sight 
distance at the cluster redshift. Angular distance is con- 
verted to proper transverse distance on the top axis. The 
dashed horizontal lines indicate the velocity dispersion of 
component A (Eq. ||), the dotted lines denote limits for 
the samples used in Sect. 2.2. The solid curves mark the 
escape velocities for masses inside radius i? of 1, 3 and 
5 X IO^'^A/q. Arrows mark the position of the giant arc 
and the X-ray and weak shear detection limits. 
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Fig. 5. Velocity dispersion profiles for galaxies with neg- 
ative and positive velocity with respect to the mean red- 
shift of cluster component A, which was held fixed at 
ZA = 0.3946. The profiles were computed using sliding av- 
erages in bins containing 30 galaxies. The error bands give 
IfT errors determined from bootstrap resampling (10000 
realisations). The curve for positive velocities was shifted 
by 300 km s^^ for clarity, with the corresponding velocity 
values marked on the right-hand axis. The bottom panel 
shows the ratio of the dispersions for galaxies with nega- 
tive and positive velocities. 
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Fig. 7. The "observed" redshift distribution of particles 
in the collision between the simulated clusters. The left 
panel shows the initial conditions whilst the right panel 
shows the data 3 Gyrs after the collision. This is to be 
compared with Fig. ||. 
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Fig. 8. The surface mass density profiles projected along 
the merger axis before (solid curve) and after the collision 
(dashed curve) . The central surface mass density after the 



collisi on is very close to that measured by Tyson et al 
(19981) . 
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Fig. 9. Comparison of rcdshift histograms of the data 
of D99 (corrected for the misidentifications discussed in 
Paper I) and Paper I (these include the data of D99), 
constrained to the area covered by D99, clustercentric dis- 
tance < 5'. Bin widths and centres are the same as in the 
inset of Figure 2 of D99. 
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